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Abstract 

The one-dimensional Hubbard model is known to possess an extended su(2) 
symmetry and to be integrable. I introduce an integrable model with an extended 
su{n) symmetry. This model contains the usual su{2) Hubbard model and has a set 
of features that makes it the natural su{n) generalization of the Hubbard model. 
Complete integrability is shown by introducing the L-matrix and showing that the 
transfer matrix commutes with the hamiltonian. While the model is integrable 
in one dimension, it provides a generalization of the Hubbard hamiltonian in any 
dimension. 
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The two-dimensional Hubbard model |I| was first introduced as a model for the descrip- 
tion of the effects of correlation for d-electrons in transition metals. The two-dimensional 
Hubbard model was then shown to be relevant to the study of high-Tc superconductivity 
of cuprate compounds. Few exact results are known and the model is still actively 
investigated. 

In contrast to the two-dimensional model, the one-dimensional Hubbard model is 
integrable. However ever since the one- dimensional model was recognized as integrable, 
its peculiar integrable structure still stands alone outside an integrable hierarchy. In this 
letter I introduce an n-state generalized model which contains the usual su{2) model. 
This model has been forecast in 0. It is studied in one dimension where it is shown to be 
integrable. Note however that the hamiltonian density is independent of the dimension 
of the lattice and therefore generalizes the Hubbard model in any dimension. 

The one- dimensional fermionic su{2) hamiltonian was first diagonalized by means of 
the coordinate Bethe Ansatz Q, a 'proof of integrability. A Jordan- Wigner transforma- 
tion shows that the fermionic hamiltonian is equivalent to a bosonic one; both forms of 
the integrable structure were investigated in the context of the quantum inverse scatter- 
ing method |^. This algebraic framework unifies integrable one-dimensional systems 
and two-dimensional classical statistical mechanics problems 0. The transfer matrix of 
the latter models provides a compact expression for all the conserved quantities of the 
quantum models. Showing that two such matrices commute directly implies the mutual 
commutation of the conserved charges of the quantum model. Moreover, the method 
provides a powerful diagonalization procedure of all the conserved charges and is a pre- 
requisite to studying the model in the thermodynamic limit with the Thermodynamic 
Bethe Ansatz. 

I first introduce the bosonic hamiltonian and the cubic conserved charge. The two- 
dimensional 'covering' statistical model and the conserved quantities are then defined. A 
calculation shows that the hamiltonian commutes with the conserved quantities. Finally 
some comments and possible directions are outlined. 

Let i?"^ be the n x n matrix with a one at row a and column (3 and zeros otherwise. 
I define the su{n) Hubbard hamiltonian as: 

H2 = E^^m + E4+i + ^E^' (1) 

i i i 

= E E {^ErE- + x-'ErE- +{E-. E')) + UY^ip^ + ^)iP^ + ^) 

where p = J2a<n E°"^ — {n — and primed and unprimed quantities correspond to 

two commuting copies of the E matrices. The complex free parameter x is a deformation 
inherited from the XX model. I am considering periodic boundary conditions. The 
hamiltonians h and h are just su{n) XX hamiltonians 0. These terms correspond to 
particle- hopping and the coupling term is an on-site Coulomb-like interaction. For |x| = 1 
the hamiltonian is hermitian. For n = 2 and x = 1, and using Pauli matrices, the 
hamiltonian is just the integrable bosonic version of the usual Hubbard hamiltonian 0: 

= ^E(^r^^i + ^f^ii) + - + f/E^f^:^ (2) 
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It is possible to write the n = 3 and x = 1 hamiltonian in terms of su{3) Gell-Mann 
matrices: 



(3) 



i a^l, 2,3,1 



2^3 



2^3' 



(3) 



I then constructed a cubic charge which commutes with H2. Investigation of the su{3) 
case showed that two copies of the only integrable models which satisfy the Reshetikhin 
criterion can only be coupled in an integrable way through A^. Another argument 
leading to a single coupling along p will be the natural generalization of the integrability 
proof. Let 



E 



E [x^ErE^.E^:^ - x-^E^E^.E^^^) 



a<n 



+ E {x-'ErEtt,E^:, - x'ErE^-,E:^, 



a,/3<n 



uUY: [{xErE-:, - x-'ErEt:,){p, + p:+, + n-2) + {Ep ^E' p) 



The conserved cubic hamiltonian is: = h^ + h^ + h'^. A direct calculation shows that H2 
and H3 commute. The shift of the generator p appearing in the coupling terms is necessary. 
The existence of a cubic conserved quantity is a strong indication of integrability; however, 
as in the 5^(2) case there is no boost giving i/3 0. 

I now construct the transfer matrix which is the generator of the infinite set of con- 
served quantities. Consider the i?-matrix of the su{n) XX model [0]: 



R{X) 



a(A) [^"" ® -^"^] 

Q,/3<n 



+c(A) E (^''^ ® + ® ^''") 



(4) 



a<n 



where a(A) = cos(A), h = sin(A) and c(A) = 1. The functions a, b and c satisfy the 
'free-fermion' condition: + b"^ = c^. Consider also the matrix 



h h 
Io{h) = cosh(-) + sinh(-) CoC'o 



(5) 



where C = J2a<n E°"^ — E"'^ and h is so far a free parameter. The transfer matrix is the 
trace of a product of L-matrices over the auxiliary space 0: 



Loi{X) — Iq Roi R-oi Iq ? t(A) — Tro (LoAf----^oi) 
where M is the number of sites. The conserved quantities are given by 

/rfPlnr(A)\ 



P+i 



dXP 



p>0 



(6) 



(7) 



A=0 



2 



where one drops trivially commuting contributions. This gives back the quadratic and 
cubic quantities written above, provided relation ( p!4D is satisfied. The non-existence of a 
boost, noted earlier, is explained by the structure of the Lax, L, matrix: L (0) — L ^(0) is 
not proportional to the identity operator (L = PL and P is the permutation operator). 

The proof that H2 commutes with r follows the lines of the first reference in |^ . This 
is one more indication that we are considering the natural multistate generalization of the 
Hubbard model. The commutator can be written as 

[H2, r(A)] = Ttq {LoM...[Hii+i, Loi+iLoi]...Loi) (8) 

i 

Let Roii+i = [Hii^i, Loi+iLoi], where Ha+i is the hamiltonian density. Let 

A tedious but straightforward calculation yields the following two equations which are at 
the root of the integrability proof: 

[hii+i, Roi+i]Roi + Roi+i[hii^i, Roi] = Roi+iDoi — Doi+iRoi (10) 
[hii+i, Roi+i]CoRoi + RQi+iCo[hii^i, Roi] = — i?oi+iCo-Doj + -Doi+iCoi?oj (H) 

Using these equations one finds 

Roii+l = -^Oi+lQoi " Qoi+lLoi , Qoi = ^oDqiIq — [CjCj,Loi] (12) 

o 

Writing the operator Q as a sum of a symmetric piece and an antisymmetric piece one 
obtains 

(13) 

The exact expression of A is not needed while that of B is similar to that of 0. A 
calculation then shows that for 



n'^U 2ab n^U 

smh.['2h) = 

one has 



sinh(2/i) = ^~^ = sin(2A) (14) 



( + 2fc2\ 

i?o.= (^^^j[^o.,/o1 (15) 

Combining this with eq. (p!3D, substituting into eq. (H), and using the cyclic structure of 
the transfer matrix we find [i/2, t(A)] = 0. 

We have shown that all the quantities obtained from the transfer matrix are conserved 
thus establishing the integrability of the model. The conserved quantities will also turn 
out to commute among themselves. One just needs to show that two transfer matrices 
at different spectral parameters commute. This follows trivially if one finds an i?-matrix 
intertwining the monodromy matrices. Then one can start diagonalizing the hamiltonians 
by the method of the algebraic Bethe Ansatz. 
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I conclude with some remarks. The eigenstate \n) of i?"" seems to play a special role. 
However the Hubbard model can be defined with respect to any other eigenstate |ao), 
and is unitarily related to the foregoing model by a unitary transformation built out of 
the transformation exchanging |t2,) and |q!o)- 

In one dimension the local symmetry is (sM(n — 1) ©'u(l)) x [[su{n — l)®u{l)). This is 
readily seen by showing that the operators X^i E^'^^ and E'^^ ^ a , (3 < n, commute with 
the transfer matrix. The primed copies are also symmetries. The proof outlined in 
holds here. Finally a transformation of the Jordan- Wigner type (when n > 2) to fermionic 
variables cannot exist because of dimensional considerations. Whether the above local 
symmetry can be extended to a non-local one involving i?"", i?"", is an open question. 

Acknowledgement: I thank P. Mathieu for fruitful discussions. 
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